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Abstract 

We show that two Pu nucleoli resonance poles near the ttA threshold, obtained in several 
analyses, are stable against large variations of parameters within a dynamical coupled-channels 
analysis based on meson-exchange mechanisms. By also performing an analysis based on a model 
with a bare nucleon state, we find that this two-pole structure is insensitive to the analytic structure 
of the amplitude in the region below irN threshold. Our results are M po i e = (1363lg — i79jj) 
MeV and (1373^Jq — i 114^g 4 ) MeV. We also demonstrate that the number of poles in the 1.5 GeV 
< W < 2 GeV region could be more than one, depending on how the structure of the single-energy 
solution of SAID is fitted. For three-pole solutions, our best estimated result of a pole near iV(1710) 
listed by Particle Data Group is (1829^gjp — i 192^^) MeV which is close to the results of several 
previous analyses. Our results indicate the need of more accurate ttN reaction data in the W > 
1.6 GeV region for high precision resonance extractions. 



PACS numbers: 14.20.Gk, 13.75.Gx, 13.60.Le 



I. INTRODUCTION 



An important task in hadron physics is to extract nucleon resonances from 7rN reaction 
data. The extracted resonance parameters are needed to understand the spectrum and 
structure of excited nucleons within QCD. They are also the starting point for analyzing 
electromagnetic meson production reaction data which have been of high precision and 
extensive in recent years 

There exists several approaches j2^14] to extract nucleon resonances (N*) from ttN reac- 
tion data. In general, almost all 4-stars nucleon resonances listed by Particle Data Group [l5[ 
(PDG) are found in all approaches. However the existence of some N* states, in particular 
those in the higher mass region, is controversial. The most investigated case is the number 
of resonances in ttN Pn partial wave. In the region near Roper A r (1440), two poles close to 



the 7rA threshold were found in Refs. [3H5|, [12] and in our recent extraction [14j, while only 



one pole in the similar energy region was reported in Refs. 0, 0, Hc| . In the higher mass 
region, the iV(1710) in PiTiN par tial- wave is not reported in Refs. [5l. Il2|. but is identified 
in all other analysis @, HMIl 

To make progress, it is important to address a commonly asked question on the extent to 
which the extracted resonance parameters depend on the reaction models employed and the 
accuracy of the empirical partial-wave amplitudes used in the analysis. For Pn resonances, 
this was investigated by Cutkosky and Wang [3| and more recently by Ceci, Svarc, and 
Zauner [9| within the Carnegie- Mellon University-Berkeley model 0|i (CMB). In an analysis 
including tcN, r)N and pseudo irnN channels, it was demonstrated [9j that the existence of 
iV(1710) depends on the structure of the nN amplitude which is related to the coupled- 
channels effects due to rjN channel. In this work, we carry out a similar investigation within 
a dynamical coupled-channels model [l6[ (EBAC-DCC). The main difference between our 
approach and CMB model is to define the non-resonant amplitude by using the meson- 
exchange mechanisms. We thus have provided additional information for examining the 
dependence of the Pn resonances on the reaction models employed in the analysis. 

Our investigation has two parts. First we would like to examine the stability of the two- 
pole structure of Pn resonances near the 7rA threshold (W ~ 1.3 GeV). Our objective is to 
examine how this two-pole structure is sensitive to the parameters of the meson-exchange 
mechanisms within the EBAC-DCC model used in our extraction }14T |. In the fits [TtJ 
(JLMS) of ttN data these parameters were determined within the ranges known from 
previous studies of meson-exchange mechanisms. Here we allow them to vary much more 
freely such that several models with different analytic properties are obtained for examining 
whether the resulting pole positions are stable within the EBAC-DCC model. 

The two-pole structure is also reported by Doring et al. [12] in an analysis based on a 

19] , the analytic 
is rather different from the 



As discussed in Ref 
13-122 



meson-exchange model with a bare nucleon state [18 
structure of this model as well as other similar models 
EBAC-DCC model, in particular in the region near the nucleon pole, mainly because of the 
differences in deriving 23( three-dimensional scattering equations from relativistic quantum 
field theory. To further examine the stability of the two-pole structure of Pn resonances and 
the existence of A(1710) state within the meson-exchange models, we also perform fit s by 
using such a model. Our formulation is similar to that developed by Pearce and Afnan 20 . 

We will show that the positions of two poles near the 7rA threshold extracted from all 
of the meson-exchange models constructed here are rather stable. This explains why the 
similar two-pole structure is also found in the other analyses jl, 0, [l2[ which use very different 
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reaction models. 

The second part of our investigation is to examine the extent to which the structure of 
the P\\ amplitude in higher invariant mass (W) region can influence the two-pole structure 
near the 7rA threshold. Here we also follow Ref. jsj to examine how the number of resonance 
states in the region near iV(1710) state listed by PDG depend on the structure of the data. 
We thus consider both the energy- dependent and single-energy solutions (SP06) of SAID jBf, 
hereafter we call them SAID-EDS and SAID-SES, respectively, as well as a solution from 
CMB [3] collaboration. The CMB amplitudes could be outdated, but is used here only for 
investigating the dependence of the Pn poles on the accuracy of the data. We will show 
that the number of resonance poles in the 1.6 GeV < W < 2 GeV could be more than 
one, depending on how the structure of the amplitude is fitted. Our results indicate the 
importance of improving the accuracy of empirical partial-wave amplitudes. More accurate 
ttN reaction data from the new hadron facilities, such as the Japan Proton Accelerator 
Research Complex (J- PARC), are needed. Our conclusion is consistent with the finding of 
Ref. 0] in which the importance of also fitting the nN — > rjN amplitude is demonstrated in 
identifying P u A(1710) state. 

In Sec. [Til we give a brief description of the coupled-channels models used in this work. 
The results are given and discussed in Sec. IHIl Section [TV] is devoted to the discussions on 
possible further developments. 



II. DYNAMICAL COUPLED-CHANNELS MODELS 



In this section we first recall briefly the EBAC-DCC model 16[ used in this work. We 
then describe how the model can be modified to obtain a model with a bare nucleon, which 
has the main feature of other ttN reaction models with a bare nucleon 12, 



A. EBAC-DCC model 

The EBAC-DCC model describes meson-baryon reactions involving the following chan- 
nels: irN, t]N, and ttttN which has ttA, pN, and aN resonant components. The ex- 
citation of the internal structure of a baryon (B) by a meson (M) to a bare N* state 
is modeled by a vertex interaction Tmb^n*- The meson-baryon (MB) states can in- 
teract via interactions v MBjM i B i which describe the meson-exchange mechanisms deduced 
from phenomenological Lagrangians. Within the model, the partial-wave amplitude of the 
M(k) + B(—k) —7- M'(k') + B'(-k') reaction can be cast into the following form (suppressing 
the angular momentum and isospin indices): 

T M B,M'B'(k,k',E) = t M B,M'B>(k,k',E) +tf IBM , B ,(k,k',E), (1) 

where the first term is defined by a set of coupled-channels integral equations 

t M B,M'B'(k,k',E) = v M B,M'B'(k,k') 

+ J2 q 2 dqv M B,M"B"(k,q)G M "B"(q,E)tM"B",M'B'(q,k',E).(2) 

M" B" M" B lt 
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Here Cmb is the integration contour in the complex-g plane used for the channel MB. The 
term associated with the bare N* states in Eq. ([!]) is 



(k',E), (3) 



where the dressed vertex function T N *^ M > B >(k, E) is calculated [17| from the bare vertex 
r ' N*^M'B'(k) and convolutions over the amplitudes t M B,M'B'(k, k 1 \E). The inverse of the 
propagator of dressed N* states in Eq. Q is 

[D-\E)] hJ = {E- m Q Nt )8 id - E id (E), (4) 
where m° N * is the bare mass of the i-th N* state, and the N* self-energy is defined by 

i 

%i,j(E) = ^2 q 2 dqr N *^ MB (q,E)G M B(q,E)r MB ^ N *(q,E). (5) 

MB M B 

Defining E a (k) = [m 2 a + k 2 ] 1 / 2 with m a being the mass of particle a, the meson-baryon 
propagators in the above equations are: Gmb^, E) = 1/[E — E M {k) — E B {k) + ie] for the 
stable iiN and r/N channels, and Gmb^, E) = 1/[E — E M {k) — E B {k) — T> MB (k, E)] for the 
unstable ttA, pN, and <jN channels. The self energy E^f^(fc, E) is calculated from a vertex 
function defining the decay of the considered unstable particle in the presence of a spectator 
re or N with momentum k. For example, we have for the 7rA state, 

y (b F\ = mA [ 2 r M wN (q) |/A^jv(g)| 2 r „x 

^rAl*. E ^ {k) q a <l [M 2 M + P] l/2 E _ En{k) _ [M 2 N{g) + k 2y/2 + 

where M n N(q) = E^q) + E]y(q) and fA->nN(q) defines the decay of the A — > ttN in the 
rest frame of A, C3 is the corresponding integration contour in the complex-g plane. The 
self-energies for the pN and aN channels are similar. 

To search for resonance poles, the contours Cmb and C 3 must be chosen appropriately 
to solve Eqs. ([2])-® for E on the various possible sheets of the Riemann surface. The 
procedures for performing this numerical task have been discussed in Ref. (13 . 24] |. Like all 
previous works (B], 10 1, we only look for poles which are close to the physical region and 
have effects on the ttN scattering observables. All of these poles are on the unphysical sheet 
of the iiN channel, but could be on either unphysical (u) or physical (p) sheets of other 
channels considered in this analysis. We will indicate the sheets where the identified poles 
are located by (s^tv, s vN , s n7! -N, s^a, s pN , s^n), where s M b and s^n can be u or p. 



B. Model with a bare nucleon state 

To examine further the model dependence of resonance extractions, it is useful to also 



perform analysis using models with a bare nucleon, as developed in, for example, Refs. [12 



20M22| . Within the formulation given in Sec. Ill such a model can be obtained by adding 
a bare nucleon (N ) state with mass m° N and iVo — > MB vertices and removing the direct 
MB — )• N — > M'B' in the meson-baryon interactions vmb,m>b'- All numerical procedures 
for this model are identical to that used in the JLMS analysis, except that the resulting 
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amplitude must satisfy the nucleon pole condition. Here we follow the procedure of Afnan 
and Pearce (20| . 

For simplicity, we include one bare Nq state and only one bare N* state. The amplitude 
can still be written in the form of Eq. ([T]) and the propagator D(E) of the term t R of 
Eq. 02]) is a 2 x 2 matrix. The nucleon pole condition can be most transparently defined by 
introducing an orthogonal matrix U T U = 1 (UJa = Uji) to diagonalize D~ 1 (E) of Eq. fll]). 
The term t^ N nN of Eq. fl3]) can then be cast into the following diagonal form 



with 



t%^ N (k,k,E)= JF^M^, 

i=l,2 

where m\ = m° N and m' = m%» , and the mass shifts are 

t^E) = i{m°-m? + £ n (£) + £ 22 (£) 

-{{ml + E 22 (E) - m? - E n (E)) 2 + A^E)} 1 ' 2 } , , 
t 2 (E) = i{m°-m° + E n (£) + E 22 (£) 

+[(m° + £ 22 (£) -ml- ^n(E)) 2 + A^E)] 1 / 2 }. 

Here Hij(E) are defined by Eq. (jSJ). The transformed vertices in Eq. (jHJ) are 

j 



(7) 



(8) 



(9) 
(10) 

(11) 

(12) 



where N* = N or N*, and the transformation operator U is defined by 

1 



22 



-C7i 



1 + ^2)1/2* 
V 



21 



with 



(1 + 1/2)1/2' 

En — Ei s 22 — s 2 



El2 E12 

Suppose E = mN pole is found in the first term of Eq. ©, we then expand 



E-m£-£i(E) = E -ml - |Ei(tojv) + 
= {E-m N )\l 



7^-£i(£) 



(E - m N ) + 



E=m* 



-S-Zi(E) 
dE u ; 



+ 



E=mpj 



(13) 
(14) 

(15) 



(16) 
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where we have defined the nucleoli pole 



Tn jy = + Ei(miv) • 



(17) 



This is the first nucleon pole condition taken into account in constructing the bare nucleon 
model. 

Defining the renormalized vertex as 



F nNN {k) = hMk)Z~ 1/2 



(18) 



with 



we then have 



tf(k -> k on , k -> k on , E -> m N ) 



[-p7rJVJv(feon)]' 
E — 777 TV 



(19) 



(20) 



Here the on-shell momentum is defined by E = ^jm 2 N + k^ n + ym 2 + k^ n . Below E = 
m/v + fn n , k on becomes positive or negative imaginary. Here we take the positive imaginary 
since we look for the physical nucleon pole. The second nucleon pole condition then defines 
the renormalized vertex F n ^^(k on ) as the physical ttNN form factor. Following the partial- 
wave decomposition procedure given in Ref. [l6j], we find 



F„ NN (k on ) = F%$j(k 



with 



fr, 



NN 



(27r) 3 / 2 77V 



V12nk, 



EN(k) + m N 



E N (k) + 777 AT 



(21) 



(22) 



where f n NN = y/4w x 0.08. Following the previous approach, the bare Nq — > vertex 
T NotnN (k) is parameterized as Eq. (I2"2"j) except that f^NN is replaced by a bare coupling 
constant f® NN , and the form factor is introduced. Explicitly, it is written as 



It 



NN 



Vl2nk A 



EN(k) + 777 /V 



1 + 



E„(k) 



(271)3/2 777^""'^ 2E N (k) ^JJ^[ E N (k)+m N 
where we use the following form factor, 

F(k,A 7rNN ) 



A: 



NN 



, k 2 + A-Inn , 



F(k,A^ NN ), 
(23) 

(24) 



The cutoff parameter A^/vat of the form factor and the bare coupling constant f® NN are varied 
along with other parameters of the model to fit the empirical 7rA^ scattering amplitudes and 
the pole conditions (fTTj) and (j5TJ. 
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TABLE I: The resonance pole positions Mr for P\\ [listed as (ReM/j, — IhiMr) in the unit 
of MeV] extracted from various parameter sets. The location of the pole is specified by, e.g., 
(snNi s v n, Sn-n-N, s n A, s p n, s a ]y) = (upuupp), where p and u denote the physical and unphysical 
sheets for a given reaction channel, respectively, xid 1S defined by Eq. (|26p . 



Model 


upuupp 


upuppp 


uuuupp 


uuuuup 




SAID-EDS 


(1359, 81) 


(1388, 83) 






2.94 


JLMS 


(1357, 76) 


(1364, 105) 




(1820, 248) 


3.55 


liV*-3p-H 


(1357, 74) 


(1363, 111) 




(1792, 280) 


2.41 


liV*-3p-L 


(1359, 69) 


(1371, 112) 




(1940, 242) 


5.33 


2iV*-3p 


(1368, 82) 


(1375, 110) 




(1810, 82) 


3.28 


27V*-4p 


(1372, 80) 


(1385, 114) 


(1636, 67) 


(1960, 215) 


3.36 


27V*-4p-CMB 


(1379, 89) 


(1386, 109) 


(1613, 42) 


(1913, 324) 


4.91 


!AT liV*-3p 


(1363, 81) 


(1377, 128) 




(1764, 137) 


2.51 



Here we note that the pole condition ( TT7l) depends on both m? and as can be seen in 
Eq. OH]) for Si(£ l ). Thus the mass renormalization of the physical nucleon includes not only 
the meson cloud effects, but also the contribution from the bare N* state. If we drop the 
iV* state, the nucleon pole condition become the usual form 



m N = 171° + Sii(mjv). 



(25) 



We use the exact conditions ( fTTl) and ( 12~TT) in our investigations. Our approach is not 
completely consistent with the rigorous approach of Ref . 2(| , but is sufficient for our present 
limited purpose of investigating model dependence of resonance extractions. Qualitatively, 
this model contains the main feature of the coupled-channels model developed in Ref. [12J 
in handling the nN scattering in P n partial wave. The main difference is in the derivation 
of meson-baryon potential vmb,m'b' from phenomenological Lagrangians, as discussed in 
Ref. [H . 



III. RESULT 

We first discuss the parameters of the coupled- channels models described in section II, 
which are varied in performing x 2 fits to empirical P\\ amplitudes using MINUIT. The non- 
resonant amplitude t MBjM / B / of Eq. Q is determined by the coupling constants and cutoffs 
of form factors of the meson-exchange interactions vmb,m'B' through solving the coupled- 
channels integral equation ([2]). In the JLMS fit [17| to the tiN data, these parameters were 
constrained within the ranges known from previous studies of meson-exchange mechanisms, 



as discussed in Ref. [16|. Here we allow them to vary much more freely such that several 
models are obtained for examining whether the resulting pole positions are stable against 
the variation of the analytic properties of the resulting amplitudes. 

In the absence of theoretical input, our main challenge is to determine the bare N* mass 
m° N * and the N* — > MB vertex function. For Pn partial wave, the number of N* parameters 
is : N N *+N N * x V" n VtMB , where N N * is the number of the bare N* and n VtMB is the number 

MB 

of parameters needed to parameterize each N* — > MB vertex function T^*^mb- 111 our fit 
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we have N N * = 1 or 2, and n vM B = 2 (MB = irN,r]N,nA,aN) or 4 (MB = pN) from 
the coupling constants guB and cutoffs Amb (as explained in Ref. We have total five 

channels (N^b = 5). We thus face a many-parameters problem in fitting the data, which is 
also present in using the CMB models with N MB = 8, 6, 3 in Refs. jjjj], [13], [§], respectively. 
We also note that the similar many-parameters problem is also a concern in all approaches 
of resonance extraction which require high precision fits of 7rN data. This common problem 
poses difficulties in assigning the errors for the determined model parameters. We thus 
follow all previous works and will only assign errors in the determined Pn resonance pole 
parameters which are determined non-linearly by the model parameters associated with 
meson-exchange interactions vmb,m'B' and bare N* states. 

Our fitting procedure is as follows. We first adjust the parameters associated with the 
meson-exchange interaction vmb,m'b' to fit Pn amplitude at low energies W ^ 1.2 GeV. To 
control the number of parameters associated with bare N* states, we then include only one 
bare N* state and try to fit the data in the entire considered energy region by adjusting 
its bare mass m° N , and vertex function parameters qmb and Amb- If this fails, we then 
also allow the parameters associated vmb,m'B' to vary. If this fails again, we then include 
one more bare N* state and repeat the process. In the region below W = 2 GeV, we find 
that the considered Pn amplitudes can be fitted with one or two bare N* states. Most 
of the resulting cutoff parameters are in the range of [5 00 MeV - 1500 MeV], which are 



similar to those in typical meson exchange models [18|, I20H22I. |25( . The bare N* masses are 
searched within the range m° N * < 2500 MeV. The interpretations of these resulting bare N* 
parameters with hadron structure calculations remain to be developed. At the present time, 
they should be considered purely phenomenologically and only the extracted resonance pole 
parameters have well defined physical meaning. 

Once a fit is obtained, we then apply the method of analytic continuation of Ref. j24j 
to find resonance poles, as also briefly described at the end of Sec. Ill Al The errors of the 
resonance parameters are then estimated by using all values obtained in all fits we have 
performed. 

For each of fits to be presented below, we assess its quality by evaluating its \ 2 P er data 
point defined by 



2 



Re[T model (Wi)] - Re[T data (iyi)] Im[T modcl (W;)] - Im[T data (Wi)] 



^ i= w^data 1 |Re(5[T d ^(^)])| 2 llm^T^iy,)])! 2 



2 



(26) 



where T data (Wj) and £[T data (Wj)] are the values and errors of the considered data, respec- 
tively; Nw is the number of the energy points where the data exist; iVdata = 2Nw is the 
number of the data points (note that there are real and imaginary components at each 
energy point). We use the single energy solution SAID-SES as data in our fits, except in 
one fit using CMB data (see below). As a reference we take the energy dependent solution 
SAID-EDS as T mode \E) to get x 2 pd = 2 - 94 , as listed in the first row of Table I along with 
their values of Pn resonance pole positions. Note that their sheet assignments are different 
from ours since they do not have aN channel. Also they do not have a pole at higher energy 
region. 



We now proceed to present our results by first recalling the three Pn poles extracted [14 
from using the JLMS parameters. They are listed in the second row of Table |T] and the 
corresponding amplitudes (solid curves) are compared with the SAID-EDS ^ (open circles) 
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FIG. 1: (Color online) The real (left) and imaginary (right) parts of the on-shell P\\ amplitudes 
as a function of the ttN invariant mass W (MeV). The solid curves are from the JLMS fit; the 
dashed (dotted) curves are from the lA*-3p-H (lA*-3p-L) fit to the SAID-EDS 0j up to W = 2 
GeV (W = 1.6 GeV); the open circles are the SAID-EDS 0]. T is unitless in the convention of 
Ref. fl. 



in Fig. [TJ Here we note that the xf,d from the JLMS fit listed in Table [I] is comparable to 
that of SAID-EDS. In general, we find it is rather difficult to get a fit with Xpd — 2-5 within 
meson-exchange model, mainly because the errors of SAID-SES are very small in W < 1.45 
GeV region within which the reproduction of the rapid sign changes of empirical amplitude 
is rather difficult due to the need of delicate balance between the attractive and repulsive 
effects in different energy region. 

In the following subsections, we present results from various fits by varying the dynamical 
content of the EBAC-DCC model as described above and using a model with a bare nucleon 
described in Sec. Ill Bl 



A. 1A*-3 P -H and lA*-3p-L fits 

We first consider the simplest variation of the JLMS fit by including only one bare N* 
state, instead of two, to fit the SAID-SES solution. In these fits, the parameters of meson- 
baryon interactions vmb.m'b' of Eq. ([2]) are taken from JLMS. We also examine how the 
extracted resonance poles depend on the data included in the fits. Here we present results 
from two fits. The solution lN*-3p-H fits the SAID-SES up to 2 GeV, while the lV*-3p-L 
to only 1.6 GeV. These two fits are compared with the JLMS results in Fig. [TJ The resulting 
resonance poles are listed in the third and fourth rows of Table [B We see that the first two 
poles near the irA threshold (~ 1.3 GeV) are in good agreement with those from JLMS. 
This suggests that these two poles are only sensitive to the data below about 1.5 GeV. 
The differences between these two fits and JLMS at higher W > 1.5 GeV mainly affect the 
positions of their third poles, as seen in Table [H 

The results presented here also indicate that with the nonresonant amplitudes of JLMS, 
only one bare N* state is sufficient to describe the ttN scattering data up to 2 GeV. All of 
the fits to be presented below are obtained by starting with non-resonant amplitudes which 
are chosen to be different from that of JLMS by tuning the parameters of vmb,ai'B'- It turns 
out that in these fits, using the procedure described above, two bare N* states are needed 
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to get comparable y 2 . 

B. 2iV*-3p and 2iV*-4p fits 

Here we investigate the dependence of the extracted resonances on the accuracy of the 
employed partial-wave amplitudes by considering the SAID-SES solution which show some 
oscillating structure in the high W ^ 1.5 GeV region. Such a structure is absent in the 
SAID-EDS (open circles in Fig. [1]). From the empirical point of view, it raises the question 
on whether the fits to the smooth SAID-EDS miss some resonance physics of the original 
ttN data. Before more precise empirical amplitudes are available, it is necessary to explore 
the extent to which these experimental uncertainties can affect the resonance extractions. 
We explore this issue by allowing the parameters associated with meson-baryon interaction 
vmb,m'B' to deviate from the JLMS values in varying these parameters along with the 
bare N* parameters in minimizing Xpd- ^ n general, the resulting 7riVA and pNN coupling 
constants from these new fits are weaker than the JLMS values and hence give rather different 
non-resonant amplitudes t n N,nN- 

We have obtained several fits which differ from each other mainly in how the oscillating 
structure of the data at high W are fitted. The results from the 2iV*-3p (dotted curves) and 
2iV*-4p (dashed curves) fits are compared with the JLMS fit (solid curves) in Fig. [2J The 
resulting resonance poles are listed in the 5th and 6th rows of Table [B Here we see again 
the first two poles near the irA threshold from both fits agree well with the JLMS fit. This 
seems to further support the conjecture that these two poles are mainly sensitive to the data 
below W ~ 1.5 GeV where the SAID-SES has rather small errors. However, the 2iV*-4p 
fit has one more pole at Mr = 1630 — i45 MeV. This is perhaps related to its oscillating 
structure near W ~ 1.6 GeV (dashed curves), as shown in the Figs. [2](b) and [2](d). On 
the other hand, this resonance pole could be fictitious since the fit 2iV*-3p (dotted curve) 
with only three poles are equally acceptable within the fluctuating experimental errors. Our 
result suggests that it is important to have more accurate data in the high W region for a 
high precision resonance extraction. 

C. 2iV*-4p-CMB fit 

To further explore the dependence of the resonance poles on the data, we consider a 
solution from CMB collaboration j3j. This solution differs significantly from the SAID-SES 
mainly at W > 1.55 GeV. For our present purpose of investigating the stability of the lowest 
two poles near the 7rA threshold, we fit the data which is obtained from replacing SAID-SES 
in the high W > 1.55 GeV region by the CMB solution. The results (dashed curves) from 
this fit with all parameters allowed to vary within the EBAC-DCC model are compared 
with that of the 2iV*-4p (solid curves) in Fig. |3j We see that both have oscillating behavior 
near W ~ 1.6 GeV and this could be the common reason why both have an addition pole 
near W ~ 1.6 GeV, as seen in rows 6 and 7 of Table [B The large differences in their fits 
at high W make their poles near W ~ 1.9 GeV very different; in particular their imaginary 
parts. On the other hand, their lowest two poles near the 7rA threshold are close to other 
fits discussed so far. This again supports the above observation that these two poles are 
determined only by the data below W < 1.5 GeV which are reproduced very well in all fits. 
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FIG. 2: (Color online) The real (upper panels) and imaginary (lower panels) parts of the P\\ am- 
plitudes as a function of the ttN invariant mass W (MeV). The JLMS (solid) results are compared 
with the results from the 2iV*-3p (dotted) and 2iV*-4p (dashed) fits. The points with errors are 
from the SAID-SES 0. T is unitless in the convention of Ref . [5(] . 



D. lA liV*-3p 

Here we consider the question concerning whether the analytic structure of the employed 
reaction model in the W < + m 7T unphysical region can strongly influence the resonance 
extractions. We first note that most of the resonances listed by Particle Data Group [l5[ 
(PDG) are from analysis which treat nucleon as a structureless basic degree of freedom 
in describing the ttN reactions; such models are used in SAID 0| and CMB [§}. Similar 
simplification is used in formulating the EBAC-DCC model [16J. On the other hand, a 
more elaborated approach has been taken to analyze irN data using models within which 
the nucleon is made of a bare nucleon A"o and meson clouds. Such models 12J, I2CH22 



need to account for the nucleon pole condition, as described in Sec. Ill B\ in fitting the 
ttN reaction data. While all of these models give similar P\\ amplitudes from threshold 
W t h = rnN + m-K to about 1.6 GeV, their analytic structure as function of the com plex energy 
could be very different in the W < + region where all dynamical models 0, E3, 0- 



22l . l25j have various singularities due to t he p arameterization of the considered meson-baryon 



interactions. This is discussed in Ref. [19(. The question is whether such differences can 
lead to very different resonance poles. 

We investigate this issue by comparing the results presented above with that from the fits 
using the model with a bare Ao described in Sec. Ill Bt In these fits, the parameters of the 
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FIG. 3: (Color online) The real (left) and imaginary (right) parts of the P\\ amplitudes as a 
function of the vriV invariant mass W (MeV). The 2A*-4p fit (solid) and the 2A*-4p-CMB fit 
(dashed) are compared with the data. The open circles with errors are from the SAID-SES [jj, 
and the crosses with errors are from the SAID-SES at W < 1.55 GeV and the CMB solution at 
W > 1.55 GeV. T is unitless in the convention of Ref. Isjl. 



meson-baryon interaction vmb,m'B' are adjusted along with the parameters associated with 
N and N* in fitting the SAID-SES up to W = 2 GeV under the nucleon pole condition (ITT)) 
and fl2~TT) . The results from one of the fits (dashed curves) are compared with the JLMS fits 
(solid curves) in Fig. HJ We see that the two fits agree very well below W — 1.5 GeV, while 
their differences are significant in the high W region, as seen in the right-hand-side panels 
of Fig. |H The resulting resonance poles are given in the last row of Table HI Similar to all 
of the cases discussed above, we also see here that the first two poles near the irA threshold 
are close to those of JLMS. Our results seem to indicate that these two poles are rather 
insensitive to the analytic structure of the amplitude in the region below ttN threshold, and 
are mainly determined by the data in the region + rn n < W < 1.6 GeV. The third pole 
from this fit is close to that of JLMS, except that its imaginary part is smaller as seen in 
the first and last rows of Table [TJ 



E. Averaged values of the extracted P\\ resonances 

To get the averaged values of the extracted Pn resonance poles, we take the values listed 
in Table I except those from models liV*-3p-L and 2iV*-4p-CMB which are obtained from 
fitting different sets of data, as described above. We further omit the values from 2A^*-4p 
in the evaluation since it has one more pole due to its oscillating behavior (dashed curves in 
Fig. 2) which needs further investigations although it is within the experimental uncertainties. 
Our values are listed in Table [Til The errors are assigned by the differences between the 
largest and smallest values listed in Table I. 

The model parameters from our fits are not relevant to the discussions given above and 
are therefore not presented. These information is available upon requests. 
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FIG. 4: (Color online) The real (upper panels) and imaginary (lower panels) parts of the P\\ 
amplitudes as a function of the ttN invariant mass W (MeV). The JLMS fit (solid) and the 
liVoliV*-3p fit (dashed) are compared with the SAID-SES (Hi- T is unitless in the convention of 
Ref. 



TABLE II: Averaged values of the extracted P\\ resonances [listed as CReMft, — ImM^) in the unit 



of MeV]. Here we identify these poles with the states listed by PDG lfij|. The location of the pole 
is specified by, e.g., (s n i^, s v n, s^v, s^Ai s pN, So-n) = (upuupp), where p and u denote the physical 
and unphysical sheets for a given reaction channel, respectively. 



States 


Location 


Averaged values (MeV) 


AT(1440) 


{upuupp) 
(upuppp) 


(1363ig,79ijj) 
(13731^ H4+14) 


iv(mo) 


(uuuuup) 


(18291^,192^0) 



IV. SUMMARY AND DISCUSSIONS 

In this work we have investigated the extraction of Pn nucleon resonances. By performing 
extensive fits to SAID-SES, we show that two resonance poles near the 7rA threshold are 
stable against large variations of parameters of meson-exchange mechanisms within EBAC- 
DCC model [16|. This two-pole structure is also obtained in an analysis based on a model 
with a bare nucleon state. Our results indicate that the extraction of Pn resonances is 
insensitive to the analytic structure of the amplitude in the region below irN threshold. 
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By performing different fits to the structure of SAID-SES as well as the old, perhaps also 
outdated, CMB data, we demonstrated that the number of poles in the 1.5 GeV < W < 
2 GeV region could be more than one. Thus our determination of the resonance poles in 
this higher W region is not so conclusive. We can only report one pole near iV(1710) state 
listed by PDG, in agreement with several previous analyses. Our results indicate the need 
of more accurate ttN reaction data in the W > 1.5 GeV region for high precision resonance 
extractions. In particular, accurate inelastic amplitudes for rjN, 7rA, pN, and aN channels 
are highly desirable for our 5-channels analysis. This will allow simultaneous fits to both 
elastic and inelastic amplitudes to firmly determine the nucleon resonances in the 1.5 GeV 
< W < 2 GeV region. The importance of performing multi-channel fits was demonstrated 
in a recent 3-channels CMB analysis [9] in which it was shown that a simultaneous fit to 
both ttN 7iN and nN — > rjN is needed to establish iV(1710) state. Thus it is important 
to obtain more extensive data of 7rN reactions including polarization observables such that 
high precision partial-wave amplitude analysis of ttN — > inrN data can be performed. Such 
experiments are possible in new hadron facility J-PARC in Japan. 

Finally, we like to mention that the analysis of electromagnetic 7r and 2n production 
data can help confirm the nucleon resonances extracted from ttN reaction data, although 
its main objective is to extract electromagnetic properties of nucleon resonances. On the 
other hand, some resonances, which have small branching ratios to it and 2tt channels and 
have large ones for KY and uN channels, could be identified from analyzing the data of 
7iV — > KY, uN which have been accumulated extensively in recent years. This is also an 
important task in N* study before the hadronic data for these channels become extensive 
at new hadron facility. 
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